Two-channel pseudogap Kondo and Anderson models: 
Quantum phase transitions and non-Fermi liquids 
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We discuss the two-channel Kondo problem with a pseudogap density of states, p(u>) oc |aj| r , of 
the bath fermions. Combining both analytical and numerical renormalization group techniques, we 
characterize the impurity phases and quantum phase transitions of the relevant Kondo and Anderson 
models. The line of stable points, corresponding to the overscreened non-Fermi liquid behavior of 
the metallic r — case, is replaced by a stable particle-hole symmetric intermediate-coupling fixed 
point for < r < r max ~ 0.23. For r > r max , this non-Fermi liquid phase disappears, and instead a 
critical fixed point with an emergent spin-channel symmetry appears, controlling the quantum phase 
transition between two phases with stable spin and channel moments, respectively. We propose low- 
energy field theories to describe the quantum phase transitions, all being formulated in fermionic 
variables. We employ epsilon expansion techniques to calculate critical properties near the critical 
dimensions r = and r = 1, the latter being potentially relevant for two-channel Kondo impurities 
in neutral graphene. We find the analytical results to be in excellent agreement with those obtained 
from applying Wilson's numerical renormalization group technique. 



I. INTRODUCTION 

The two-channel Kondo effect represents a prime ex- 
ample of non-Fermi liquid behavior arising from a sta- 
ble intermediate-coupling fixed points Theoretically, its 
physics is essentially understood, thanks to an exact solu- 
tion by Bethe ansatz^i 3 - In addition, boundary conformal 
field theory (CFT) has proved to be a powerful technique 
to study the low-energy properties of the multi-channel 
Kondo model^ allowing, in particular, the calculation 
of exact asymptotic Green's functions^ CFT techniques 
have also been used to calculate exact crossover Green's 
functions^ Further, by means of Abelian bosonization 
and subsequent rc-fermionization it has been possible to 
map the two-channel Kondo problem onto a resonant- 
level model which reduces to a free fermion form for a 
particular value of exchange anisotropy.— 

On the experimental side, a number of heavy-fermion 
materials, displaying deviations from Fermi-liquid behav- 
ior, have been speculated to realize two-channel Kondo 
physics arising from non-Kramers doublet ground states 
of U or Pr ions4~— However, to our knowledge, there is 
no unambiguous verification of these proposals. Con- 
sequently, various attempts were made to realize the 
two-channel Kondo effect in nanostructures, and indeed 
success was reported^ for a setup of a semiconductor 
quantum dot coupled to two reservoirs Very recently, 
signatures of two-channel Kondo behavior of magnetic 
adatoms on graphene have been reported^ and this mo- 
tivates to discuss two-channel Kondo impurities in non- 
metallic hosts. In particular, neutral graphene realizes 
a pseudogap density of states (DOS), p(uS) cx \uj\ r with 
r = 1, at low energies. 

The single-channel pseudogap Kondo problem has 
been studied extensively in the context of Kondo impu- 
rities in unconventional superconductors. The main dif- 



ference to the familiar metallic Kondo problem^ is the 
absence of screening at small Kondo coupling J, leading 
to a quantum phase transition upon increasing JJ£rlH 
The universality class of this phase transition changes as 
function of and the relevant low-energy field theories 
have been worked out in detail in Refs. I20il2ll . 

Although two-channel Kondo physics has been specu- 
lated about in the context of graphene^ the two-channel 
pseudogap Kondo model has received little attention. A 
central question is about the fate and character of the 
non- Fermi-liquid phase at finite r. To our knowledge, 
the only study of the model has been reported in a brief 
section of Ref. [l9|, but there only numerical results were 
given for small bath exponents r£2- 

The purpose of this paper is to close this gap: We shall 
investigate the two-channel Kondo and Anderson models 
with a pseudogap DOS in some detail, using both ana- 
lytical and numerical renormalization group (RG) tech- 
niques. Our main findings for the two-channel Kondo 
model are: 

(A) The overscreened non-Fermi liquid (NFL) phase 
of the metallic two-channel Kondo model^ survives for 
< r < r max ~ 0.23, albeit with an important mod- 
ification: It is no longer represented by a line of NFL 
fixed points (where particle-hole (p-h) asymmetry is 
marginal), but instead there is only an isolated stable 
p-h symmetric NFL fixed point, i.e., p-h asymmetry is 
irrelevant for r > O^ 3 . Furthermore, this stable NFL 
phase is only reached for Kondo couplings larger than 
a critical coupling, i.e., a boundary quantum phase tran- 
sition emerges between a local-moment phase with an un- 
screened spin moment and the NFL phase. In contrast, 
for r > r max the NFL fixed point disappears, leaving 
only two stable phases with unscreened spin or channel 
(i.e. flavor) moment, respectively, which are separated by 
a quantum phase transition. 
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(B) The two-channel pseudogap Kondo physics for 
both r < 1 and r > 1 can be fully understood in the lan- 
guage of the two-channel Anderson model, by virtue of a 
generalization of the approach presented in Ref. [U The 
low-energy field theory describing the quantum phase 
transition between the phases with free spin and flavor 
moments is given by a level crossing of a spin doublet and 
a flavor doublet minimally coupled to conduction elec- 
trons. Similar to the single-channel pseudogap Kondo 
problem, r = 1 is found to play the role of an upper- 
critical dimension, where the hybridization is marginal. 
For r > 1, the transition is a level crossing with perturba- 
tive corrections, whereas a non-trivial critical fixed point 
emerges for r < 1. This fixed point is shown to display 
an emergent spin-channel Z2 symmetry. As in the single- 
channel case, none of the quantum phase transitions is 
described by a Landau-Ginzburg- Wilson-type theory of 
a bosonic order parameter, instead all are "fermionic" in 
nature. 

The following subsection gives a more detailed sum- 
mary of our results. 

A. Summary of results 

The two-channel Kondo model with a pseudogap host 
density of states can be written as H. = Hk + %t>, with 

n K = ^ [jK5 ? -cL(o)wc^(o) + y cL(o)c lCT (o)}i) 

i 

,A 

T-Lb = \^ / dk \k\ r kc\ ia Ckia- 
I J 

Here, we have represented the bath, Hb, by linearly dis- 
persing chiral fermions Ckia, where i = 1, 2 is the channel 
index. S is a spin-1/2 SU(2) spin, f is the vector of Pauli 
matrices, summation over repeated spin indices a is im- 
plied, and Cj CT (0) = J dk\k\ r Ckia is the conduction elec- 
tron operator at the impurity site. The spectral density 
of the Ci a (0) fermions follows the power law \uj\ r below 
the ultra-violet (UV) cutoff A; details of the density of 
states at high energies are irrelevant for the discussion 
in this paper. In addition to the Kondo coupling Jk, we 
have also included a potential scatterer of strength Vq at 
the impurity site which will be used to tune p-h asym- 
metry. (An asymmetry of the high-energy part of the 
DOS would have a net effect similar to non-zero Vo; for 
simplicity we will assume in the following that the DOS 
is p-h symmetric.) 

As we shall show below, a comprehensive analysis re- 
quires to consider - in addition to the two-channel Kondo 
model - the two-channel Anderson model, commonly 
written as H = Ha + T~Lh with£ 

(J i 

+ 90 E(i ff ) (*i c «* + H - c -) + V oJ2 4(°)^(°)- 

kai ai 



Here, the isolated impurity has four states, i.e., a spin 
doublet I cr) = I f), I I) and a channel doublet \i) = 
|1), |2). Their mass difference, e = £s — £q, will play 
a role as a tuning parameter of the quantum phase tran- 
sition. 

For a given value of the bath exponent r, the two- 
channel pseudogap Kondo and Anderson models display 
common RG fixed points. The phase diagram and crit- 
ical behavior depend on r, with r = 0, r = r max , and 
r = 1 marking qualitative changes and playing the role 
of critical "dimensions" . In the following, we describe 
our central results for the phase diagrams and RG flows, 
which are partially consistent^ with the ones reported 
in Ref. [nj. The qualitative behavior is visualized in the 
RG flow diagrams in Fig. [1] for the two-channel Kondo 
model and Fig. [5] for the two-channel Anderson model, 
respectively. In the latter case, a cut through the RG 
flow at Vo = is shown. 

The metallic case r = 0, has been studied extensively, 
and a line of infrared stable NFL fixed points governs the 
behavior at any finite coupling - this is the well-known 
two-channel (or overscreened) Kondo effect. In the two- 
channel Anderson model, this line of fixed points can 
be accessed by varying So, i.e., initial parameters with 
different £0 flow to different fixed points along this line^ 5 . 
- note that this flow leaves the v = plane for e ^ 
(dashed in Fig. [2j all symbols denote the renormalized 
coupling parameter). 

For positive r with < r < r max , the line of stable NFL 
fixed points collapses to an isolated p-h symmetric NFL 
fixed point. In addition, the local- moment fixed point 
(LM) of an unscreened spin moment now becomes stable. 
In the language of the Kondo model, LM corresponds to 
j = v = 0, while in the Anderson model it corresponds 
to s = —00, g = 0. The phase transition between LM 
and NFL is controlled by a critical p-h symmetric fixed 
point (SCR); note that this p-h symmetric fixed point is 
located outside the v = plane for the Anderson model 
shown in Fig. As r — > SCR approaches LM and 
the critical behavior of SCR is pcrturbatively accessible 
for small Kondo coupling Jk . The phase diagram of the 
Anderson model is mirror symmetric, i.e., there exists 
also an unscreened channel (or flavor) local-moment fixed 
point LM' at e = 00 , g = 0^1 and a corresponding critical 
intermediate-coupling fixed point SCR' at positive e. 

As r — > ?*max SCR approaches NFL, and the two fixed 
points disappear for r > r max . In the p-h symmetric 
Kondo model, this implies that the flow is towards LM 
for any value of j, but for large asymmetries, LM' may 
be reached^ In the Anderson model, the hybridization 
g remains relevant at e = for r max < r < 1, but the 
flow is towards a single unstable intermediate coupling 
fixed point (ACR) in the v = plane, i.e., ACR is p-h 
asymmetric, but is invariant under the Z2 transformation 
Q. At finite coupling, the transition between the two 
stable fixed points LM and LM' is controlled by ACR. 

Finally, as r — > 1 ACR moves toward g — > and for 
r > 1 the phase transition becomes a level crossing with 
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FIG. 1: Schematic RG flow diagrams for the two-channel 
Kondo model with a pseudogap DOS, p(co) oc \uj\ r ■ The hor- 
izontal axis denotes the renormalized Kondo coupling j, the 
vertical axis is the renormalized potential scattering v, repre- 
senting particle-hole asymmetry. Dashed flow lines symbolize 
a flow out of the plane shown here. The thick lines correspond 
to continuous boundary phase transitions; the full (open) cir- 
cles are stable (unstable) fixed points, a) r — 0, i.e., the famil- 
iar metallic case. For any finite j the flow is towards the line 
of NFL fixed points, describing two-channel non-Fermi liquid 
behavior, b) 0<r<r max : P-h asymmetry is irrelevant in the 
NFL phase, such that the line of fixed points is replaced by 
a single p-h symmetric NFL fixed point. The local-moment 
fixed point LM is stable and separated from NFL by a p-h 
symmetric critical SCR fixed point. For r — > (r — ^ r max ), 
SCR approaches LM (NFL). Depending on microscopic de- 
tails, a second NFL' phase may be reached at large couplings 
and asymmetries, separated by a critical ACR fixed point, see 
text, c) r>r max : The NFL phase disappears, and the only 
phase transition is between LM and LM', the latter repre- 
senting a free channel (i.e. flavor) moment. This transition is 
controlled by ACR. Note that the character of this transition 
changes at r — 1, where ACR merges with FImp, see text. 



perturbative corrections, controlled by the free impurity 
fixed point (FImp) at e = 0, g = 0. 




FIG. 2: As Fig. [T] but for the two-channel Anderson model. 
The horizontal axis denotes the energy difference e between 
spin and flavor impurity levels, the vertical axis is the renor- 
malized hybridization g. The diagrams represent cuts, taken 
at v = 0, through the full RG flow. The flow diagrams are 
mirror-symmetric by virtue of the Z2 transformation (Q. a) 
r — 0. The line of NFL fixed points, describing two-channel 
non- Fermi liquid behavior crosses the v = plane at e — 0. 
b) 0<r-<r- max : The fixed points LM, LM' are stable, corre- 
sponding to unscreened spin and flavor moments, respectively. 
Within the v = plane, NFL is replaced by ACR, whereas 
two isolated p-h symmetric NFL/NFL' fixed points exist out- 
side this plane. The transitions to NFL/NFL' are controlled 
by the p-h symmetric SCR/SCR' fixed points (located out- 
side the v = plane), c) r max < r < 1: g is still relevant at 
s = 0. However, NFL is now replaced by a single unstable 
fixed point (ACR) located in the v = plane. At finite g, the 
transition between the two stable fixed points LM and LM' 
is controlled by ACR. d) r > 1: g is irrelevant, and the only 
transition is a level crossing (with perturbative corrections) 
occurring at g = s = 0, i.e., at the free-impurity fixed point 
(FImp). 



B. Outline 

The bulk of this paper is organized as follows: We start 
in Sec. |TT] by discussing the relevant impurity models, 
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suitably generalized to higher degeneracies, along with 
their underlying symmetries. In Sec. Mil we present se- 
lected results from Wilson's numerical rcnormalization 
group (NRG) for the two-channel pseudogap Anderson 
and Kondo models which illustrate the content of the 
flow diagrams in Figs. Q] and [2] In particular, we show 
properties of the non-trivial intermediate-coupling fixed 
points as function of the bath exponent r. Sees. II VI and IVl 
are devoted to the epsilon expansion studies of the criti- 
cal fixed points, using the variables of the Kondo model 
(Sec.|lV| and that of the Anderson model (SecEJ. The 
latter provide access to the physics near the upper-critical 
dimension r = 1. Concluding remarks will close the pa- 
per. A discussion of the Majorana representation of the 
two-channel Kondo problem and its fate in the presence 
of a pseudogap DOS is relegated to the appendix. 



II. MODELS, SYMMETRIES, AND MAPPINGS 

A. Anderson model 

The two-channel Anderson model can be understood 
as describing the level crossing between two impurity 
doublets - one spin doublet and one channel (i.e. flavor) 
doublet - coupled to conduction electrons via a hybridiza- 
tion term. The model features an SU(2) spin x SU(2)fl avor 
symmetry. This can be straightforwardly generalized to 
SU(A r ) sp in x SU(if )flavor symmetry, where N is the num- 
ber of spin degrees of freedom and K the number of fla- 
vors. The Hamiltonian can be written as: 

H A = U h +e s ^2\a)(a\+e q J2\a){a\ (3) 

a a 

+ goY,(W)(a\c ka a +H.c.) + V&^cL, (0)^(0) 

kaa oca 

Here, the conduction electrons c kaa transform under a 
fundamental representation of SU(AT) and SU(iT) and 
carry the corresponding spin a and flavor a indices, a in- 
dicates a transformation behavior according to the conju- 
gate representation. For N = 2, K = 1, the Hamiltonian 
in Eq. §3§ describes the single-channel Anderson model in 
the limit of infinite Coulomb repulsion (U = oo), studied 
using RG in Ref. [2l|. In case of a metallic host, r = 0, 
the multi-channel Anderson model is integrable and has 
been solved using the Bethe Ansata^r— and the numer- 
ical RG method^ 

The metallic two-channel Anderson model, i.e. N = 2, 
K = 2, has been proposed as a model for the ob- 
served non-Fermi liquid behavior of the heavy-fcrmion 
superconductor UBei3^ In this scenario, the 51 2 ground 
state of the U ion is identified as the F3 non-magnetic 
quadrupolar doublet, while the first excited state is the 
5f 3 magnetic doublet.— This then can promote a 
quadrupolar Kondo effect where the quadrupolar dou- 
blet is quenched by the hybridization with conduction 
electrons which carry both magnetic and quadrupolar de- 



grees of freedom^ In particular, since the energy differ- 
ence between the two doublets appears to be small, a 
mixed valence state is likely requiring the study of the full 
Anderson models Consequently, the model with a 
pseudogap DOS is of potential relevance not only to two- 
channel impurities in graphene, but also to quadrupolar 
Kondo impurities in unconventional superconductors. 

The Anderson model §3§ is not particle-hole symmet- 
ric for any value of Vb, due to the asymmetric structure 
of the impurity. However, p-h symmetry is dynamically 
restored for < r < r max both inside the NFL/NFL' 
phases and at the critical SCR/SCR' fixed points, see 
Fig.S 

Interestingly, for N = K and a p-h symmetric bath, 
the Anderson model displays a spin-channel symmetry, 
i.e., is invariant under the combined transformation 

\a) O \a) 
e O -£o 

V o -Vb (4) 

Here, the spin-carrying impurity states are transformed 
into the flavor-carrying states and vice versa, i.e., the 
two SV(N) sectors are interchanged, together with a p-h 
transformation. 



B. Kondo models 

The Anderson model has two Kondo limits. On the 
one hand, for Eq — £ s — e q — > — oo it maps to a _ftT-channcl 
SU(A f ) S pi n -Kondo model, where a spinful impurity is cou- 
pled to K channels of conduction electrons. For N = 2 
the Hamiltonian reads 

n K = n h + j k s ■ cL(o)w<w(o) (5) 

+ Vo5>L(0)c« CT (0) 

oca 

where S is a spin-1/2 SU(2) spin and a, a' =t, \.. For > 
2 the impurity spin is in a fundamental representation 
of SU(A). The parameters of the Kondo model (O are 
related to that of the Vq = Anderson model ^ through: 

Jk = NV = & . (6) 

Nl 

The Kondo limit is reached by taking Eq — > — oo, go — > oo, 
keeping Jk fixed. Note that a potential scattering term 
is always generated. 

On the other hand, for e +°° the Anderson 
model can be mapped to a A^-channel SU(A')fl avor -Kondo 
model, where 5* represents a SU(A') impurity which is 
screened by the N spin degrees of freedom of the conduc- 
tion electrons. Such multi-channel flavor Kondo effect is 
relevant, e.g., to the charging process of a quantum box, 
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where the flavor degree of freedom is taken by the phys- 
ical chargei 31 ' 32 

We note that the multi-channel Kondo model cannot 
be obtained by a Schricffer- Wolff transformation from 
any standard Anderson model (i.e., written with free- 
electron operators and local Coulomb interaction). 



C. Large- TV limit 

The SU(TV) multi-channel Kondo model can be solved 
in a dynamic large-TV limit for both fully symmetric 
(bosonic) and fully antisymmetric (fermionic) representa- 
tions of SU(A r )^ 3 - The fermionic version of this solution, 
with K = 7 TV and K,N — > oo, has been generalized to 
the pseudogap case^ The large- TV phase diagram, Fig. 1 
of Ref. [HI, is similar to that of the N = K = 2 case, 
i.e., the overscreened non- Fermi liquid phase survives for 
small r, where it is reached for a certain range of cou- 
plings only, while this phase disappears for larger r < 1. 
Also, all leading anomalous dimensions vanish for r > 1. 

Two qualitative differences between the large-TV sce- 
nario and TV = K = 2 arc worth noting: (i) The critical 
"dimension" r max of TV = K = 2 splits into two in the 
large- TV limit, with their values and the detailed behav- 
ior depending upon the value of 7. (ii) The quantum 
phase transitions in the large-TV limit are governed by 
lines of fixed points, with continuously varying exponents 
as function of the particlc-holc asymmetry, in contrast to 
the isolated critical fixed points SCR and ACR in Figs.Q] 
and[2j This implies that the scaling dimension of Vq van- 
ishes at criticality as TV — »• 00, rendering the large- TV limit 
partially singular. Wc therefore refrain from a detailed 
discussion of the models © and (|5|) for large TV. 



D. Observables 

In the bulk of the paper, we will focus on a few im- 
portant observables which characterize the phases and 
phase transitions of the impurity models under consid- 
eration. Those include the correlation- length exponent, 
the impurity entropy, various susceptibilities, and the 
conduction-electron T matrix (or impurity spectral func- 
tion). Their definition is standard, and we refer the 
reader to Rcfs. Il9ll2lll35|-[37l for a detailed exposition. 
Here we only summarize a few key aspects. 

Spin susceptibilities, x^ spm \ are obtained by coupling 
external magnetic fields both to the bulk and impurity 
degrees of freedom as explained in detail in Ref. [33- For 
the impurity part, here, this reads 

-H im rtW)KA°'\, (7) 

where -ffi mp! i is the magnetic field coupling to the impu- 
rity spin, while the \ l aa , with i — 1,TV — 1 are genera- 
tors of SU(TV). In the following, we exploit the SU(TV) 



symmetry and only evaluate the corresponding suscep- 
tibility tensor in the 1-direction choosing the represen- 
tation A 1 a , = g(5 CT ,i<5 CT ',i — ^cr, 2#<t', 2)-— We proceed as 
usual by calculating the magnetic susceptibilities via the 
corresponding linear response functions. Note that the 
impurity susceptibility is composed of 

Ximp(T) = Ximpamp + 2Xb,imp + (Xb,b - Xb,b)> (8) 

where Ximp.imp is the response to H lmp , Xb,b measures the 
bulk response to the field applied to the bulk, Xb,imp are 
the cross terms, and Xb b denotes the bulk response in the 
absence of the impurity. Flavor susceptibilities, x^ flavor \ 
can be defined in the Anderson model in analogy to the 
spin susceptibilities (i.e., with a — > a). 

Owing to symmetries, the total magnetization in both 
the spin and flavor sectors is conserved. This implies that 
the impurity contributions to the spin and flavor suscep- 
tibilities, and Ximp V ° r ' ) : do n °t acquire anomalous 
exponents at the intermediate-coupling fixed points, but 
instead obey Curie laws with (in general) fractional pref- 
actors. In contrast, the local spin and flavor susceptibil- 
ities follow anomalous power laws, x^ 1 "' 1 <x T~ 1+r, * P 

and Xioc™^ x T _1+, 'x , with universal r-dependent 
exponents r\ x . We note that a direct calculation of both 
susceptibilities is only possible in the Anderson model, 
as the Kondo limit suppresses the local piece of one of 
the susceptibilities. To shorten notation, we employ the 
convention x = x^ spm ^ and rj x = r)^ in the following. 

Similar to T'x^ > p 11 ' flavor \ the impurity entropy ap- 
proaches a universal fractional value as T — > 0. The 
conduction-electron T matrix, on the other hand, follows 
an anomalous power law similar to the local susceptibil- 
ity, T(w) cx cj- 1+7 i t . 

At the non-Fermi-liquid fixed point of the familiar 
metallic two-channel Kondo model (r = 0), power laws 
are replaced by logs, XlooXimp hil/T - this also im- 
plies that the prefactor of the leading Curie term in Ximp 
vanishes due to an exact compensation. 



III. SELECTED NUMERICAL RESULTS 

The NRG technique 3 ^ is ideally suited to study prop- 
erties of quantum impurity models, including non-Fcrmi 
liquid phases and quantum phase transitions. Initial 
NRG results for the two-channel pseudogap Kondo model 
were shown in Ref. Il9l . Here we extend and comple- 
ment this early analysis by NRG results for the two- 
channel (TV = K = 2) Anderson model. We per- 
form explicit calculations for a bath density of states 
p{uj) = (1 + r )/(2D)\uj/D\ r e(D 2 - lo 2 ) with D = 1. 
Unless otherwise noted, we employ NRG parameters 3 ^ 
A = 6 and TV, = 600. 

The qualitative behavior of the two-channel Anderson 
model is summarized in the flow diagram in Fig. [2J In 
addition to the stable LM/LM' phases, these flow di- 
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FIG. 3: NRG for the impurity susceptibility, Tximp, at the 
intermediate-coupling fixed points ACR (□), SCR (A), and 
NFL (0)- Also shown are the results from the renormalized 
perturbation theory of the Kondo model, which allows to ac- 
cess SCR near r = [Sec. HV] Eq. Jl4]), dashed], and that of 
the Anderson model, appropriate for ACR for r < 1 [Sec. [V] 
Eq. (|49[1 , dash-dot] . The SCR data have been partially taken 
from Ref. 0. The vertical dashed lines indicate the critical 
dimensions r = r max and r = 1. 



agrams feature three non-trivial fixed points: the sta- 
ble NFL/NFL' fixed points and the critical fixed points 
SCR/SCR' and ACR. Some of their key properties are 
summarized in Figs. [3] and |4l which show the numeri- 
cally determined impurity contributions to the spin sus- 
ceptibility and the entropy, respectively, together with 
analytical results obtained from the epsilon expansion 
of Sees. [TV] and El These plots nicely show that NFL 
and SCR approach each other as r — > r max , while ACR 
evolves continuously near r max . The fixed-point proper- 
ties also show that SCR approaches LM as r — > 0, with 
^Ximp — > 1/4 and Si mp — > In 2, and ACR approaches 



1, with Tx 



imp 



-> 1/8 and Si 



FImp as r 

Further, the stable NFL fixed point follows Tx 



imp 
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FIG. 4: As Fig. O but for the impurity entropy 5i mp . The 
perturbative expansions are in Eq. (|15p (dashed) and Eq. (|53p 
(dash-dot). The dotted horizontal lines correspond to Si mp = 
0.5 In 2, In 2, and In 4. The inset shows a zoom onto the region 



FIG. 5: Phase diagram of the two-channel Anderson model 
as function of the energy difference eo and the bath exponent 
r, keeping <?o = 4 fixed. The discontinuous change at r = r max 
is apparent. 



and S'imp — >0.51n2asr— > - the well-known properties 
of the metallic two-channel Kondo problem. 

The disappearance of both NFL and SCR upon in- 
creasing r beyond r max implies a discontinuous evolution 
of the phase diagram as function of r. In Fig.[5]wc present 
a cut through the phase diagram of the Anderson model 
at fixed go which illustrates this fact. We note that such 
a discontinuous evolution occurs if a stable intermediate- 
coupling fixed point disappears (here NFL); in contrast, 
if a trivial fixed point changes its nature from stable to 
unstable, the evolution is continuous, like in the single- 
channel pseudogap Kondo model. 

The correlation-length exponents v are displayed in 
Fig. [6l illustrating that r ~ and ?' max play the role 
of lower-critical dimensions for the p-h symmetric tran- 
sition controlled by SCR, with v — > oo, whereas r = 1 is 
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FIG. 6: Inverse correlation-length exponent 1/u obtained 
from NRG for the ACR (□) and SCR (A) critical points, 
together with the analytical RG results from the expansions 
in r [Sec.[lVl Eq. O, dashed] and in (1-r) [Sec.|V] Eq. flU), 
dash-dot]. Note that v of ACR becomes very large for small 
r - this corresponds to the extremely slow flow from ACR to 
NFL for < r < r max . 
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FIG. 7: NRG results for a) the impurity susceptibility Tximp 
and b) the impurity entropy Simp of the two-channel Ander- 
son model at r — 0.22. Parameter values are go = 2 and 
£o = (dashed), 0.5 (solid), 1 (long-dash), 2 (dash-dot). Both 
£o = 0.5 and 1 are located in the NFL phase while eq = 2 cor- 
responds to the LM phase. The so = 0.5 curves show the 
flow from ACR to NFL upon lowering T, where in particu- 
lar the impurity entropy increases along the RG flow. (The 
flow is too slow to reach its fixed point within the accessible 
temperature range.) 



the upper-critical dimension of the ACR transition, with 
v = 1 for all r > 1 (and logarithmic corrections at r = 1). 

Fig. [7] shows the temperature evolution of both Tximp 
and S- lmp for r = 0.22, i.e., slightly below r max . Here, the 
flow from ACR to NFL (compare Fig.[]J>) is nicely visible 
at eo = 0-5 where both Tximp and 5i mp increase along the 
RG flow. (The large value of v at ACR renders the flow 
very slow.) Remarkably, the fact Sacr < <5nfl violates 
so-called g-theorem^S which states that the impurity en- 
tropy should decrease along the flow. As this theorem 
applies to conformally invariant systems only, we con- 
clude that the fixed points under consideration are not 
described by a conformally invariant theory. (The same 
conclusion can be drawn for the quantum phase tran- 
sitions of the single-channel pseudogap Kondo problem, 
see Ref-lU Also, such "uphill flow" may occur in models 
with long-range interactions, see e.g. Ref. l4ll .) 



IV. WEAK-COUPLING EXPANSION FOR THE 
MULTI-CHANNEL KONDO MODEL 

In this section we review the standard weak-coupling 
expansion^ for the multi-channel SU(2)-Kondo model 
in Eq. ([S]), extended to a pseudogapped bath density of 



states! 15 ' 19 As we show below, this expansion captures 
the properties of the critical fixed point SCR at small 
r. In principle, it also allows to access the stable NFL 
fixed point, but this requires a particular limit of large K 
which does not allow to extract quantitative results for 
the K = 2 case of interest. 

RG equations can be derived, e.g., using the field- 
theoretic schem o 43 ' 44 where logarithmic divergencies, oc- 
curring for r = 0, are replaced by poles in r by means of 
"dimensional" regularization. Doing so, r will only enter 
in the bare scaling dimension of the couplings. To two- 
loop order, the equations for the renormalized couplings 
j and v read 

- <•./ r • f ./ ''• (9) 

j3(v) = rv, 

where K is the number of equivalent screening channels. 
Importantly, there is no rcnormalization of v, a result 
which persists to higher orders. 

Apart from the LM fixed point, j = v = 0, the function 
in Eq. ((9]) yields two further zeros, given by j = 
(I ± VT — 2Kr)/K. The smaller one corresponds to an 
infrared unstable fixed point at 

j* = r + yr 2 +C(r 3 ), v* =0 (10) 

which can be perturbativcly controlled for r — > and 
any K. We label this fixed point by SCR. As r -> 
SCR approaches LM. The larger zero predicts an infrared 
stable fixed point at 

f = l-r- I ±r 2 + <D(r 3 ), v* =0. (11) 

Strictly speaking, this fixed point is perturbatively acces- 
sible only if the limits r — > and K — > oo are either taken 
in this order (this corresponds to r = 0) or together such 
that Kr is kept fixed. For r = 0, where v is marginal, 
this zero of j3(j) is commonly associated with the line of 
stable non-Fermi liquid fixed points of the multi-channel 
Kondo model, which exists for all K > 2. For r > now 
v becomes irrelevant, in agreement with our numerical 
results^ 3 - which show that the NFL line of fixed points 
shrinks to a single p-h symmetric NFL fixed point, FigJI] 

A. Observables near criticality 

The properties of the p-h symmetric critical fixed point 
SCR, existing for < r < r max , can be determined in 
analogy to the single-channel case, with explicit calcula- 
tions given e.g. in Ref. l37l Expanding the beta function 
© around the fixed point value (fit)]) yields the correla- 
tion length exponent v. 

i=r-^r 2 + G(r 3 ). (12) 
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The leading-order pcrturbative corrections to the im- 
purity susceptibility and entropy are given by 



A(T X 



AS, 



\ - -El 

np) — . j 

3tt 2 In 2 



imp 



Kf r. 



(13) 



The corresponding values of Tx^™^ follow via Eq. (|4]) 
from LM-H-LM'. The hybridization term, g , is irrele- 
vant at LM and LM' for r > 0. 



B. Hybridization expansion and upper critical 
dimension 



Inserting the fixed-point value j* (|10)) we obtain 
Tximp = \(l-Kr)+0{r 2 ), 



(14) 



3tt 2 



ln2 1+— Kr 6 )+0(r b ). (15) 



The anomalous exponent of the local susceptibility is 
given by rj x = Kj to leading order, which evaluates 
to 



7] x = Kr 2 + C(r 3 ) . 



(16) 



Finally, the T matrix exponent is t\t = 1 — j (with no 
factor of K, as the T matrix describes the scattering of 
electrons from one specific channel), resulting in 



T)t = 1 — r. 
Note that this result is exacts 



(17) 



V. HYBRIDIZATION EXPANSION FOR THE 
MULTI-CHANNEL ANDERSON MODEL 

We now turn our attention to the multi-channel Ander- 
son model. As we show below, the variables of the Ander- 
son model will allow us to obtain an essentially complete 
understanding of the multi-channel pseudogap Kondo ef- 
fect both for r < 1 and r > 1. A similar conclusion was 
reached for the single-channel case in Refs. I20II21T and - 
on a technical level - our work represents a generaliza- 
tion of the calculation for the infinite-C/ Anderson model 
in Ref. [H 



A. Trivial fixed points 

For vanishing hybridization go, the multi-channel An- 
derson model features three trivial fixed points: for 
£0 < the ground state is the spinful TV-fold degenerate 
local-moment state (LM) and, analogously, for £0 > it 
is the A"-fold degenerate flavor local-moment state (LM'). 
In these cases the impurity entropy equals In N and In K, 
respectively. For £0 = there arc (N+K) degenerate im- 
purity states, we refer to this as the free-impurity fixed 
point (FImp), with entropy ln(A^ + K). The impurity 
spin susceptibilities are 



(spin) _ 



2F, LM 
FImp 



imp 



2(N+K) 

LM' 



(18) 



In the following we perform an expansion around the 
FImp fixed point, i.e., around £0 = 0, go = 0. 

The impurity states are represented by bosonic op- 
erators b\ for a = 1,...,K and fermionic operators 
f a for (7 = 1,...,N. Single occupancy of the local- 
ized levels is enforced by the Hilbcrt space constraint 
Q = J2 a b&ba + ^2 a flfa = 1 which will be implemented 
using a chemical potential Ao — > oo. Observables are 
then calculated a a 45 i 46 



(O) = lim 

An— >a 



(QOh 

(Q)a ( 



(19) 



where (■ • • }\ denotes the thermal expectation value in 
the presence of the chemical potential Ao- 

Furthermore, we need to introduce chemical-potential 
counter-terms which cancel the shift of the critical point 
occurring in perturbation theory upon taking the limit of 
infinite UV cutoff. Technically, this shift arises from the 
real parts of the self-energies of the b a and f a particles. 
We introduce the counter-terms as additional chemical 
potential for the auxiliary particles, 



(20) 



The SXbj have to be determined order by order in an 
expansion in go . Note that counter-term contributions in 
observables in general enter both numerator and denom- 
inator in Eq. (|T9"|) . 

In the path integral form the model © is written as 



S 



(It 



fAdr + Xo+e s + SX f )U 



b a (d T + A + £, 
go (fabaC a cr(0) + 



h 5X b )b a 
c.c.) 



dk \k\ r c kaa {d T + k)c kaa 



-A 



(21) 



where Ao is the chemical potential enforcing the con- 
straint exactly. Here, we implicitly sum over a and a. 
The model (|21[) shows a transition driven by variation 



of £ 



At the £o = go = fixed point, tree level 



scaling analysis shows that 

dim [g ] = 
dim[£ ] = 



(22) 



1 . 



As in the single-channel model ) 20 ' 21 this establishes the 
role of r — 1 as an upper-critical dimension where go is 
marginal. 
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(a) (b) 



FIG. 8: Feynman diagrams entering the self energies up 
to quadratic order in go. Full/dashed/wiggly lines denote 
fa/ccr/ba propagators, the full dots are go interaction ver- 
tices. 



diagram does not exist due to the directed nature of the 
propagators). The beta function 



dg_ 

(If i 



(31) 



can now be obtained by taking the logarithmic \x deriva- 
tive of Eq. ([25]). Since — wc can solve for /3(g) 
and finally obtain 



m _ (32) 



We perform a field-theoretical RG analysis using the 
minimal subtraction scheme . 43 ^ 44 Rcnormalizcd fields 
and dimensionless couplings are introduced according to 



to one-loop order. One can also consider the flow away 
from criticality, i.e. the flow of the renormalizcd tuning 
parameter e using S 2 insertions. The resulting correla- 
tion lengths exponent is given in Eq. ([38]) below. 



u = 




(23) 




b a = 




(24) 


C. r max < r < 1 


9o = 


V Z f Z b 


(25) 


For r < 1 the trivial fixed point g* 



where \i is the renormalization group energy scale. No 
renormalizations are needed for the bulk fermions as their 
self interaction is assumed to be irrelevant. The RG is 
conveniently performed at criticality, i.e. we assume that 
Eq is tuned to the critical line and set e q = e s = 0. 

To determine the RG beta function /3(g) we evaluate 
the fermionic self-energy up to one-loop order 



^u( iu} n) = Kg 



iu)' 



dk\k\ r - 



1 



(26) 

corresponding to the diagram in Fig. Here, we have 
introduced the abbreviated notation iuj n = iw n — Ao for 
the Matsubara frequencies iui n = iir(2n + l)/f3. In the 
limit Ao — > oo and (3 — > oo we obtain 



E/,(iw n ) - Kg 2 /i 2f [ dee r -^- 



(27) 



1 



w -Kg 2 (u?r—+ioj n ^ ) . (28) 
\ r 2r J 

An analogous expression is found for the bosonic self en- 
ergy depicted in Fig. [5Jd. The renormalization factors 
are determined such that they cancel the ^ pole in the 
self-energies minimally and render the inverse Green's 
function finite. We thus find 



1 - K 



.9 



2f' 



1 — N 



2f 



(29) 



The mass counter-terms are given by the real parts of the 
self-energies 



SX f = Kq 2 ^ — , 
r 



Ng 2 ^- 
r 



(30) 



To one-loop order, there is no vertex renormalization of 
<?, hence we have Z g =l at this order (note that a g 3 



is unstable, 

and the critical properties are instead controlled by an 
interacting fixed point (labelled ACR) at 



1 



9 



K + N 

with anomalous field dimensions 

d In Zi, 



Vb 



T)f 



Kg) 

m 



dg 
d\nZ f 



*2 



g* 



dg 



Ng 



Kg*' 



(33) 

(34) 
(35) 



The corresponding RG flow diagram is displayed in 
Fig. [2b. 

ACR describes a quantum phase transition below its 
upper-critical dimension. As a result, low-energy observ- 
ables calculated at and near ACR will be fully universal, 
i.e., cutoff- independent, and hyperscaling is fulfilled. 



D. 



r > 1 



For all r > 1, i.e., above the upper-critical dimen- 
sion, the phase transition between LM and LM' is now 
controlled by the non-interacting FImp fixed point at 
g = e = 0. Hence, for all r > 1 the phase transition 
is a level crossing with pcrturbative corrections - this re- 
sults e.g. in a jump of the order parameter Tx\ oc (see 
below), i.e., the transition is formally of first order. Con- 
sequently, hyperscaling is violated, and all obscrvablcs 
will depend upon the UV cutoff. 

For the marginal case, r = 1, we expect a logarithmic 
flow of the marginally irrelevant hybridization g, charac- 
teristic of the behavior at the upper-critical dimension. 
The RG beta function 



a, s K + N 3 



(36) 



10 



(a) 



o 



(b) 



o 



(a) 




(b) 





FIG. 9: Feynman diagrams entering (xioc)a up to quadratic 
order in go- Open circles are sources, and the blank boxes 
denotes the counter terms SX f . Notation otherwise as in Fig. [5] 



can be integrated (recall (5(g) = 
f the UV cutoff A 

? 2 W = 



where I describes 



the reduction of the UV cutoff A — s- £A) to give 

,2 



1 - [K + N)gl ln£ 



(37) 



with g(£=l) = go- This result can be used to determine 
logarithmic corrections to observables. 



E. Observables near criticality 

1. Correlation-length exponent 

We start with the correlation-length exponent, v 1 of 
the ACR fixed point. This exponent describes the van- 
ishing of the characteristic crossover temperature in the 
vicinity of the critical point T* oc (e — e*) v . The lowest- 
order result for is, which can be obtained either using the 
field-theoretic RG scheme via composite operator inser- 
tions or using the familiar momcntan shell scheme, is 



1 



v 



r + 0(f 2 ) (r<l). 



(38) 



For r > 1 the transition is a level crossing, formally 
u=l. 



2. Local susceptibility 

The local susceptibility xioc = Ximp,imp at the critical 
point follows the scaling behavior xioc T~ 1+I,x with 
an anomalous exponent r\ x . To obtain the corrections 
to the tree- level result xioc oc T" 1 we introduce a xioc 
rcnormalization factor Z x from which one obtains the 
anomalous exponent according to 



m ^ 

dg 



(39) 



We determine Z x by calculating (xioc) directly using per- 
turbative corrections up to quadratic order in g . The 
corresponding diagrams entering (xioc)a are given in 
Fig. [5J In terms of the renormalized coupling constant g 



FIG. 10: Feynman diagrams entering the corrections to the 
unperturbed part of {Q)\ up to quadratic order in go. Nota- 
tion as in the previous figures. Black boxes denote the counter 
terms <5Af,. 



we find at the energy scale /i = T 

Kg 



(Xl. 



e T 




T 



dx- 



2 tanh ■ 



9 ? 1 

X , X X 

tanh x 

2 2 2 



(40) 



Furthermore, the denominator (Q)\ receives corrections 



(Q) Ao = (N + K)e~^ 
+2NKg 2 e~^ 



2 resulting 


in 




A 












/ dx 


tanh 


--ll 


Jo x 




2 J 



(41) 



The local susceptibility xioc can then be directly obtained 
by Eq. (|19[) . The renormalization factor Z x is then de- 
termined, using minimal subtraction of poles, in an ex- 
pansion in r as 



1 - Kg 2 



1 



(42) 



and from this we can directly deduce the anomalous ex- 
ponent of the local spin suceptibility 



'•X 



2Kg* 



(fia 



2K 



: (l-r). 



(43) 



K + N 

follows by the replacement 



The expression for rj x 
N. 

Above the upper-critical dimension, r > 1, we simply 

have tj x = and, thus, , Xioc™ oc or ^ u> ~ 1 - 

For the marginal case r = 1, a calculation analogous to 
that in Ref. [2U using Eq. (|37|) . gives 



(spin) 



\\ nL0 \2K/(K+N) 



(r = l). 



(44) 



3. Order parameter 

Inside the stable phases LM and LM', xioc can be used 
to define an order parameter for the quantum phase tran- 
sition: Txioc" is finite (zero) for e < (eo > 0), similarly 
^Xb™' is finite (zero) for £ > (eg < 0). Approaching 
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the critical point, both order parameters vanish continu- 
ously according to 



(spin) 



(11a 



(X Sn 



(45) 



for r < 1, which follows e.g. from hyperscaling. Note that 
these order parameters display a jump upon crossing the 
transition for r > 1. 



4- Impurity susceptibility 



(a) 




O 



(b) O 



o 



FIG. 11: Further Feynman diagrams entering corrections to 
the impurity susceptibility to quadratic order in go- Notation 
as in the previous figures. 



The evaluation of the impurity susceptibility Xi mp to 
second order in go requires the summation of further dia- 
grams as depicted in Fig. 1111 In terms of the renormalized 
coupling g we obtain 



2(Xb,imp)A — 



Kg 2 a 
T 



-4 tanh : 



x^ 



x + 



cosh 



(46) 



and 



(Xb,b>A - (Xb,b>A - 



Kg 2 _; 
— 6 



dx- 



2 tanh ■ 



x tanh | 



cosh 



2 cosrT 



2 
(47) 



Collecting all contributions to Ximp to second order in 
g the poles present in the xioc diagrams cancel and the 
remaining momentum integrals are UV convergent for 
r < 1. Performing these integrals for r < 1 the impurity 
susceptibility reads 



Txl 



spin) 
imp 



i 



2(N+K) 
+0(g 4 ). 



1 ~g 2 K l + ln4- 



2N 
N+K 



In 4 



(48) 



As above, the expression for T\^^ follows by the re- 
placement K N. With the value of the coupling at 
the ACR fixed point fl33j) wc finally find for N = K = 2, 
to leading order in (1 — r), 



(spin) 



±(l-r)+0(f 2 ) (r<l) 
(r>l) 



(49) 



with T\ 



(fla 



j, (spir 
A imp 



due to the emergent Z2 sym- 



metry ([4]). A comparison to NRG data is in Fig.O Note 
that Tximp receives only weak additive logarithmic cor- 
rections at r = 1; multiplicative logs as in xi oc are absent 
here. The same applies to Si m p below. 



5. Impurity entropy 

The impurity contribution to the entropy can be 
obtained from the thermodynamic potential f2i mp by 
Si m p — ~ ^r^irnp- At the FImp fixed point the entropy 
is Si m p — In (TV + K), and the lowest-order correction is 
computed by expanding the thermodynamic potential in 
the renormalized hybridization g. Note that this correc- 
tion vanishes for r > 1, as g* = there. Here, we write 
the partition function in the physical sector of the Hilbcrt 
space (Q = 1) as^i 



lim 



(Qh 



Zimp.O Ao^oo (Q) Xo ' 



(50) 



where (. . -)a .o is the expectation value in the presence 
of Ao without coupling to the bath. The thermodynamic 
potential is given by 



imp,0 



-Tln Zimp 



imp,0 



(51) 



The correction to (Q) \ due to the coupling to the bath 
up to quadratic order in g has already been calculated 
in Eq. (|4T|) which enables us now to directly evaluate 
Eq. (|51[) . Taking the temperature derivative of the re- 
sulting expression and evaluating the remaining integral 
in the limit T — > and r — > 1 wc obtain 



Oii 



ln{K + N)- g 2 



ANK 
N + K 



In2 + C( 5 4 ). (52) 



As expected, the entropy correction is fully universal and 
finite in the limit T — >• 0. Note that in higher-order terms 
of the diagrammatic expansion for the thermodynamic 
potential ^i mp disconnected diagrams appear^ 

Finally, inserting the fixed-point value of the coupling 
g into Eq. (|52[) . we find the impurity entropy in the two- 
channel case N — 2 and K — 2 to be 



Si, 



ln4- (1 -r)ln2 + C(f 2 ) (r < 1) 



In 4 



(r> 1) 



(53) 



Comparison with the numerical result again shows very 
good agreement, Fig. SJ 
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6. Conduction- electron T matrix 

For the single-channel pseudogap Kondo problem, it 
has been shown that the conduction-electron T matrix 
T(uj) displays a power-law divergence of the form T(ui) cx 
|w|~ r at all intermediate-coupling fixed points3£i2i An- 
alytically, this follows - for all perturbative expansions 
- from the diagrammatic structure of the T matrix (or, 
alternatively, from a Ward identity) . 

In the two-channel case, we find the same qualitative 
arguments to apply, i.e., at the NFL, SCR, and ACR 
fixed points the T matrix obeys the exact result 

T(u) oc \oj\- r (r < 1). (54) 

At one- loop level, an explicit calculation gives tjt = (K + 
N)g 2 which yields ijt = 1 — r as expected. For r > 1, 
ImT(u) oc 8(uj). 

For r = 1, the logarithmic flow of the coupling, 
Eq. (|37p . can used to deduce T(uS) oc l/(w| lnw|) which 
gives 

ImTH0C Mln^ (55) 
in analogy to Ref. [U 



VI. CONCLUSIONS 

We have explored the two-channel Kondo effect for 
magnetic impurities embedded into a fermionic host with 
a power-law pseudogapped density of states. We have 
determined the phase diagram as function of the DOS 
exponent r and discussed the boundary quantum phase 
transitions of the relevant Kondo and Anderson models. 
These transition are described by fermionic (as opposed 
to usual bosonic) quantum field theories; from their prop- 
erties we conclude that there is no underlying CFT de- 
scription. 

Our results demonstrate the versatility of the 
Anderson-model epsilon expansions developed in 
Refs. [20II21I : Those have not only allowed a full under- 
standing of the critical behavior of the single-channel 
(S = 1/2) pseudogap Kondo problem, but also of the 
corresponding underscreened^i and overscreened pseu- 
dogap Kondo models (this work). Further applications, 
e.g., to multi-impurity models, appear possible. Also, 
the Anderson-model formulation should enable studies 
of non-equilibrium dynamics of pseudogap Kondo 
problems. 

Our results are of potential relevance to two-channel 
impurities in unconventional superconductors and in 
graphene; for the latter case the extension of the present 
calculations to finite bias^ is an interesting future topic. 
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Appendix A: Compactifled o-r Kondo model and 
0(3)-symmetric Anderson model 

Here we briefly discuss an alternative formulation of 
the two-channel Kondo model which eventually leads to 
a theory of non-interacting Majorana fermions. 



1. Metallic bath, r— 

The low-energy physics of the standard two-channel 
Kondo problem (r = 0) has been argued to be equivalent 
to that of the so-called cr-r Kondo model - this is a "com- 
pactifled" single-channel Kondo model where the roles of 
the two screening channels are taken by spin and a charge 
pscudospin . 49 i 50 The corresponding Hamiltonian can be 
expressed as 

,A 

U a - T = [Ji<?(0) + J 2 r(0)] -S + / dkkc\ a c krT (kl) 

J-A 

where, as above, S is a spin-1/2 SU(2) spin and we 
have represented the bath by linearly dispersing chiral 
fermions Cfe CT . Spin degrees of freedom a =t,i are im- 
plicitly summed. The conduction electron spin ct(0) and 
pseudospin r(0) are defined as 

3(0) = (4(0),cl(0))-f.^(0)^ (A2) 

r(0) = (4(0),c,(0)).f. (A3) 

where c a (0) = J dkeka- Interestingly, the low-energy 
physics of the cr-r Kondo model is described by a fixed 
point with non-Fermi liquid behavior which is located 
at strong coupling, not at intermediate coupling as in 
the two-channel Kondo problem. The equivalence of the 
two-channel Kondo model and the cr-r Kondo model has 
been established using bosonization and conformal field 
theory techniques 42r— 

The nature of the low-energy non-Fermi liquid be- 
comes transparent by considering the so-called 0(3)- 
symmetric Anderson model which displays an anoma- 
lous hybridization term. Its Hamiltonian is given by 
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%0(3) = UlcA + "Hahyb with 



%lcA = U [riff ~ 



(A4) 



H 



f dkk c\ a c ka + g V [/t C(7 (0) + H.c] , 

/M(°) + /! c i(°)+ H - c -l . ( A5 ) 



ahyb 



where /J creates the localized impurity state with spin 
a and nj CT = f^f a . Note that the chemical potential on 
the impurity site has been chosen such that the model is 
p-h symmetric. In the Kondo limit, the 0(3)-symmetric 
Anderson model maps onto the <t-t Kondo model.— 

Appealing to the adiabatic continuity between the 
U = and large- U limits in this Anderson model sug- 
gests to discuss the weakly-interacting case. The Hamil- 
tonian %o(3) can be conveniently re-written in terms of 
Majorana fermions: 



H 



0(3) 



" /-A 

Udi d 2 d 3 d + i ^ / dkktp- ka ipka (A6) 

3 

+igo X! i>a{0)d a + i(g - 2g a )ip (0)d . 



Here, the impurity Majorana fermions are defined by 

/t = A=(d! - ida), h = A=(~ds + ido), (A7) 

where d^ a = d a and {d a ,dp} = 6 a p for a = 0,1,2,3. 
The Majorana fermions for the conduction electrons are 
defined similarly^ which in Fourier space reads 



Cfcf 



1 

V2 



{-itpki ~ ^k2) , c kir 



1 



V2 



(A8) 



with ipl = ip-ka- Remarkably, for 2g a = go the impurity 
couples only via d a for a = 1,2,3 to the conduction 
Majorana fermions while do remains free. 

The Hamiltonian in Eq. (|A6[) is in particular suitable 
to study thermodynamic quantities. Note that for U = 
the model is exactly solvable and the impurity Green's 
functions G a {r) = — (T T d a (r)d a {Q)) are known exactly. 
For U = and 2g a = go, their Fourier counterparts read 

G {iuJn) = — , G a (iu„) = — p — r, (A9) 



for a = 1, 2, 3. Here, we have introduced Aq = ng^- 

A straightforward calculation now shows then that 
there is a residual impurity entropy S- imp = ^ In 2 of a 
free Majorana fermion^ By adiabatic continuity, this 
entropy persists into the regime of large U and then cor- 
responds to the entropy of the overscreened two-channel 
Kondo impurity. 

2. Pseudogap bath 

The obvious question is whether the o-t Kondo and 
0(3)-symmetric Anderson models continue to represent 
the physics of the two-channel Kondo problem for a pseu- 
dogap bath DOS with r > 0. To answer this, let us con- 
sider the non- interacting limit and 2g a = go of the 0(3) 
Anderson model. The Green's functions in Eq. (|A9[) for 
u> n / A <C 1 are now given by2i 



G (iu) n ) 
G a (iuj n ) 



1 



iu n + iA sgn(uj n )\uj n \ r 



(A10) 
(All) 



for a = 1, 2, 3. They yield an impurity entropy of 



1 3 
Simp = -ln2 + -rln2. 



and an impurity susceptibility of 

3r 



(A12) 



TXimp(T) 



32 



(A13) 



This result is not in agreement with the numerical data 
in Figs. [3] and SI which instead are well fitted by Si mp = 
±ln2 + 2rln2 and T Xiia p(T) = § (Ref. M). We are 
forced to conclude that the low-energy behavior of the 
cr-T Kondo and 0(3)-symmetric Anderson models is not 
identical to that of the two-channel Kondo model once 
r > 0. In other words, the equivalence is restricted to the 
metallic r = case. Given the fact that neither bosoniza- 
tion nor CFT appear to be applicable to the pseudogap 
Kondo models, this may not come as a surprise. 

Let us finish with the remark that an extension of 
the Majorana resonant-level model that corresponds to 
the solvable point of the two-channel Kondo modcli to a 
pseudogap DOS also does not yield the numerically found 
impurity entropy. 
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